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VARIATIONAL PROBLEM 
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Abstract. In this work we prove a weak Noether type theorem for a class of 
variational problems which include broken extremals. We then use this result 
to prove discrete Noether type conservation laws for certain classes of finite 
element discretisation of a model elliptic problem. In addition we study how 
well the finite element scheme satisfies the continuous conservation laws arising 
from the application of Noether's 1st Theorem (E. Noether 1918). 

We summarise extensive numerical tests, illustrating the conscrvativity of 
the discrete Noether law using the p-Laplacian as an example. 

1. Introduction 

Symmetries are an extremely important and continually occurring feature of 
differential equations arising from many applicable areas, including mathematical 
physics, meteorology and differential geometry, that was first developed by Sophus 
Lie for the purpose of studying solutions of differential equations in the late 19th 
century |Lie71[ c.f.]. 

Noether's (1st) Theorem |Noe71| is a striking result which in the continuous 
setting connects these symmetries with conservation laws associated to the Euler- 
Lagrange equations of a variational problem. Roughly, the theorem states that 
given a variational problem with an underlying symmetry, there exists a natural 
conservation law associated to it. For example, a symmetry of translation with 
respect to the spatial coordinates results in a conservation of linear momentum, a 
symmetry of rotation results in conservation of angular momentum and a symme- 
try of translation with respect to the temporal coordinate gives a conservation of 
energy. A famous example from meteorology is that of potential vorticity. This is a 
conservation law arising from a particle relabelling pseudo-group symmetry. This 
quantity is extremely important in studying the evolution of a cyclone. 

The work of Noether has gained public attention recently with the publication 
of an article in the New York Times |Angl2| where the result is 

"consider[ed] ... as important as Einstein's theory of relativity; it 
undergirds much of today's vanguard research in physics, including 
the hunt for the almighty Higgs boson." 
In the discrete setting, Noether's Theorem has been studied in terms of differ- 
ence equations |Dor011 IHM04j , where it was shown that a discrete equivalent of the 
conservation law holds when a discrete symmetry was built into the discrete La- 
grangian. In this work we turn our attention to the finite element method (FEM). 
FEMs form one of the most successful numerical methods for approximating the 
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solution to partial differential equations (PDEs) [BraOli lBS94[ ICiaTSl c.f.]. A topic 
which has been the subject of much ongoing research is that of constructing FEMs 
which inherit some property of the continuous problem. The notion of discreti- 
sations inheriting some geometric property from the continuous problem can be 
seen as a generalisation of geometric integration |HLW06[ c.f.] to the case of PDEs 
and is a rapidly developing area of research. Some of the properties studied in the 
discretisation of PDEs are the same as in the geometric integration of the ODE, 
for example the Hamiltonian structure of a given problem. Others are based on 
completely new notions, for example, the recent development of the discrete ex- 
terior calculus [AFWlOj . which, as the name suggests, is a discrete equivalent to 
the Cartan based exterior calculus. This has allowed for a rigorous description of 
discrete differential forms and the associated discrete function spaces as a discrete 
differential complex. This provides a framework which may be used as a first step 
in the construction of a variational complex in a similar light to that developed in 
[HM04j for difference equations. A first step in this direction was taken in [Man06| . 
A review of some of the huge quantity of topics arising from this area, including Lie 
group integrators, discrete gradient methods as well as FEMs for differential forms 
is given in [CMKOllJ. 

As opposed to geometric integrators, the term used for numerical methods with 
some geometric property of an ODE, the methods for PDEs arc generally called 
mimetic methods. The class of FEMs which fall under the mimetic framework are 
the mixed methods, for example the Raviart-Thomas scheme |RT77] . It is not only 
FEMs which fall under the category of mimetic methods, in fact there are finite 
difference (FD) |BBL09| and finite volume (FV) schemes which are characterised 
as mimetic. Note that also there is an intrinsic relationship between each of them. 
For example, an appropriate choice of quadrature for the Raviart-Thomas finite 
element scheme results in the mimetic FD scheme |CMR09( . 

The Lagrangian piecewise polynomial FEM is not a mimetic method. Most stan- 
dard methods cannot inherit geometric properties of the continuous PDF. There 
is an underlying algebraic condition which must be satisfied for these properties to 
be inherited by the approximation scheme HM041 IMan06) . 

The purpose of this paper is to show that variational numerical problems have 
their own conservation laws which derive from the same principle as that discovered 
by Noether, giving rise to discrete (numerical) forms of conservation laws which are 
automatically preserved by the scheme. 

The classical Noether theorem is only applicable to classical solutions of the 
variational problem. As such we derive weaker versions of the theorem applicable 
to a wider class of solutions to the problem, including the broken extremals. We will 
discuss how these laws are naturally passed down to the Lagrangian finite element 
scheme and hence quantify the discrete Noether quantity associated to this FEM. 
That is, we write the exact Noether quantity for this discretisation, in the same 
spirit as |HM04j . 

We will also study how well the Lagrangian finite element scheme satisfies the 
strong conserved quantities arising from Noether's Theorem measured in an ap- 
propriate weak norm. That is, we consider how well this finite element scheme 
approximates the Noether conservation law for the continuous problem (when one 
exists). We will also present some interesting numerical results, quantifying the 
deviation of the approximation in terms of a computable estimator which we are 
able to use to construct an h-adaptive scheme (local mesh refinement) aimed at 
minimising the violation of the smooth conservation law to a user specified toler- 
ance. 
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The paper is set out as follows: In |2]we introduce some fundamental notation 
and the model problem we consider. In ^|3]we briefly describe Noether's Theorem 
and the background material needed. To illustrate its application we apply the 
Theorem to a simple model problem. In fj4]we weaken the invariance criterion on 
which the classical Noether Theorem is based, ultimately allowing us prove two ver- 
sions of the theorem applicable to weaker solutions of the problem. In ^Jsjwe discuss 
how the results of ^ can be passed down to give discrete counterparts to our weak 
Noether's Theorem. We perform numerical experiments to demonstrate that the 
quantities derived are indeed conserved at the discrete level. We also discuss trivial 
Lie group actions (those of translation with respect to the dependent variable) and 
how the mimetic methods relate to this case. Finally, in fj6]wc study the properties 
of the finite element solution with respect to the original (strong) Noether Theo- 
rem. We also detail an interesting numerical result by constructing a computable 
estimator, aimed at measuring the violation of the strong Noether Theorem in a 
weak norm for the Lagrangian finite element scheme. We perform some numerical 
experiments demonstrating that there is a superconvergence of the estimator over 
the finite element approximation of the solution to the Euler-Lagrange equations. 
We then proceed to test an adaptive scheme based on the estimate allowing us 
to minimise the discrete violation of the continuous conserved quantity up to user 
specified tolerance. 



2. Notation 

Let J7 C R'^ be a bounded domain with boundary dfl. We begin by introducing 
the Sobolev spaces |Cia781 IEva98j 

Lp(r2) = |(/) : J |(/)|^ < ooj for p <E [1, c») and Loo(f^) = {0 : esssupQ |(/)| < cx)} , 

(2.1) 

Wp(rj) = Lp{n) : D"(/) e Lp(rj), for |a| < k} and H*^(17) := W^i^l), (2.2) 
which are equipped with the following norms and semi-norms: 

MUn)--- /j-r (2.3) 
Mi, ■■= ll-llwnn) = E l|D"«ll^,(o) (2-4) 



\a\<k 



P - u,i^ - l|D°-rL,(o) (2-5) 

\a\=k 

hWl — llf'llHno) = ll«llw,^(o) ' (2-6) 

where a — {ai, . . . , ad} is a multi-index, \a.\ — 'Y^f^iOti and derivatives D" are 
understood in a weak sense. We pay particular attention to the cases k — 1,2 and 

Wlin) := closure of C^{n) in W^(f^). (2.7) 
Let L = L(x, u, Vit) be the Lagrangian. We will let 

/[■]: ^ R (2 8) 

^ ^[0] :=^i(a;,0,V0)da;. 

be known as the action functional. The problem arising from the calculus of varia- 
tions is to seek a function extremising the action functional. For simplicity we will 
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consider the minimisation problem, that is, to to find u G Wp(r2) such that 

inf ^[v]. (2.9) 

Note that we are imphcitly coupling the minimisation problem with homogenous 
Dirichlet boundary conditions. 
We will use the notation that 

V dxi ' ' dxd / 

denotes a column vector of spatial derivatives of the Lagrangian, i.e., derivatives 
with respect to the independent variables. The derivative with respect to the de- 
pendent variable is denoted ^ and let p = (pi, . . . jPd)^ = Vu then 

dL fdL{x,u,Vu) dL{x,u,Vu)y 



9(Vm) V dpi ' ■ ■ ' ' dpd 

denotes the vector of derivatives of L with respect to the gradient of u componen- 
twise. 

It is well known [Eva98J IGH96[ c.f.] that if w is a (smooth) minimiser of the 



variational problem (2.9) then it solves the quasilinear, second order, variational 



PDE called the Euler-Lagrange equations 

2.1. Assumption (null Lagrangian). We assume that the Lagrangian is not null, 



that is, the Euler-Lagrange equations (2.121 are not satisfied by all smooth func- 
tions. 

3. Noether's First Theorem 

For the reader's benefit we will briefly describe Noether's first theorem in the 
continuous, smooth case and necessary background material. We assume, in this 



section, that L is smooth and the minimisation problem (2.91 has a solution (not 
necessarily unique) which is at least C^(r2), i.e., smooth enough to satisfy the Euler- 



Lagrange equations (2.121 



3.1. Definition (one-parameter group). The transformation 

{x, u) — 7^ (S(a;, u; e), (f>(a;, u; e)) —: {x, u) (3.1) 

is said to be a one-parameter group if the following conditions hold 

(1) The parameter choice of e = yields the identity, i.e., 

{x,u) =(S(a;,M;0),$(a;,u;0)). 

(2) The inverse is given by the parameter — e, i.e., 

[x, u) = (S(5, u; -e), $(S, u; -e)) . 

(3) The transformation is closed under composition, i.e., if 

(S, u) = (S(S, m; (5), $(S, S; 5)) 

then 

(S, u) — (S(a;, u\e + 5), $(a;, u\e + 5)) . 
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3.2. Definition (infinitesimal). The infinitesimals, ^{x,u) and (t){x,u) of the one 
parameter group are defined as 

\ V dS(a;,M;e) 

E{x,u) inn ; (3.2) 

e->0 de 

d$(x,u;e) 

(p{x,u) :— lim (3.3) 

e-s-o de 

3.3. Definition (characteristics). We define the characteristics, which are given in 
terms of the infinitesimals of the group, to be 

Q{x, u, Vu) := (j){x, u) - {^{x, u))"^ Vm. (3.4) 

3.4. Definition (variational symmetry). Let F := {{x,u{x)) : a; G T} be the graph 
of u over a subdomain such that T C ft. Also let — S(r; e), then the transfor- 
mation (3.1 ) is said to be a variational symmetry if 

holds for any smooth subdomain T of fJ. 

3.5. Theorem (infinitesimal invariance |01v93[ Thm 4.12]). A variational symme- 
try group with infinitesimals ^, (p and characteristics Q of the action functional 

^{u\ = / L(a;,u,Vu) (3.6) 

satisfies 

Proof See |01v93l Thm 4.12] □ 

3.6. Theorem (Noether's 1st Theorem j01v93[ Thm 4.29]). Suppose the variational 
problem {2.9] is invariant under the action of a one-parameter group of symmetries 
with characteristics Q. Then Q is also a characteristic of a conservation law of 
the Euler-Lagrange equation ( |2.12 1, that is, there exists a vector valued functional 
^ = ^[u\ such that 

divC^M) = Q^[u]. (3.8) 
Hence over solutions of the Euler-Lagrange equations ^\u\ = 0, we have that 

div(^[w])=0. (3.9) 

For the problem we consider in this work, that of a hrst order Lagrangian, the 
conservation law, ^ , takes the form 

dT 

nu]=Li+——M-i-^Vu). (3.10) 

Proof Using the result of Theorem |3.5| we have that 

Fit BT 

Noting by the product rule that 



then it holds that 



df ( BT \ BT 

J- dL BL\ ^. / 3L 

-'^^^a(w) + a^J+H^a(v^ + ^^. 



(3.13) 
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Since the quantity 



J- dL dL ^ 

div — — + = 
o(yu) au 



(3.14) 



are the Euler-Lagrange equations associated with the action functional (2.8) they 
vanish, conchiding the proof with 



dL 

d{Vu) 



as required. 



(3.15) 
□ 



3.7. Remark (the form of '^). It is clear from the identity (3.8) that for our model 
problem '^/f = ^/f{u,Wu). 



3.8. Remark (the beauty of the theorem). What makes Theorem 3.6 truly re- 
markable is its constructive nature. For completeness we will give an example of 
the construction of for the Laplacian. 

3.9. Example (Laplace's problem). Let us consider the case / — /(|a;|) then the 
Lagrangian, 



L{x,u,Wu) ^|Vu|' 



fu, 



(3.16) 



is invariant under the rotational group SO{d). For simplicity we restrict d = 2, we 
calculate the infinitesimals from the group of rotations, note that in this case $ = 
and 

Xi cos (e) — X2 sin (e) 
Xi sin (e) — X2 cos (e). 



E.{x,u;e) 



(3.17) 



It then holds that 



lim 



d5(a;;e) 
de 



"X2 



(3.18) 



2^2 ((c*!"")^ — (92u)^) /2 — Xidiud2U + X2fu 
Xi ({diu)^ — {d2u)j /2 + X2diud2U — xifu 



(3.19) 



(3.20) 



In this case the characteristic of the group of rotations is 

— X2diU + Xid2U. 

Making use of Theorem |3.6| we see 
is a conservation law over solutions of ^[u] = 0. 

3.10. Remark (trivial Lie group actions). For any variational problem, the Euler- 
Lagrange equations, as already mentioned, are given in variational (or divergence) 
form. As such, if we assume that L does not depend on u, that is L = L(a;, Vm), 
then the Euler-Lagrange equations themselves can be a Noether conservation law. 
Indeed, consider the case of Example 3.9 with / = 0. It is clear by definition that 
Au = div(VM) = is a conservation law. It arises from Noether's Theorem under 
the trivial Lie group action, that of translation in the dependent variable 



{x, u) — > {x, u + e) 



(3.21) 



For this action, the infinitesimals are ^ = and 0=1. 
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4. Noether's theorem for weak solutions 



Noether's Theorem (Theorem 3.6 ) as it is stated in Sj3| only makes sense for 
classical solutions of the Euler-Lagrange equations ( 2.12[ ). We wish to "weaken" 



the theorem such that it is applicable to extremals which are W^(f2) with jump 
discontinuities, the so called broken extremals jGH96| . We begin by defining a 
weaker invariance condition than that of Definition 13.41 



4.1. Definition (weak variational symmetry). The transformation (3.11 is said to 
be a weak variational symmetry if 

/ L(x,u,Vu) dx = / L(x,u,Vu] dx (4.1) 
holds over the domain fl. 

4.2. Remark (strong symmetry weak symmetry). We note that any strong 
variational symmetry is also a weak symmetry but the converse is not true. 

4.3. Tiieorem (Noether type conserved quantities for weak variational symme- 
tries). Suppose that the variational problem |2.9| ) has a weak variational symmetry 
and that the minimiser is smooth, i.e., u G C (57). Let (j) and ^ be the infinitesimal 
generators of the symmetry as in Dehnition \3.2\ Then 



. + —6 + div if) - - 



= 0. (4.2) 

4.4. Remark (structure of ( |4.2[ )). The weak conservation law given in Theorem 
|4.3| has a very clear structure. The first two terms represent the weak Euler- 
Lagrange equations, which do not vanish in general since the infinitesimal does 
not necessarily belong to the correct space of test functions. The last two terms 
essentially represent the weak conservation law itself. 

Proof of Theorem |4.3| The proof of this is based on computing the inner and 
outer variations with respect to the infinitesimals of the one parameter group (see 



Definition 3.2) |GH961 c.f.]. The inner variations are with respect to the indepen- 
dent variables and the outer variation with respect to the dependent variables. 
We define the function o : IR+ — > IR as 



o(e) J^[u + e4>]^ / L{x,u + e(t),Vu + eVcj}). 
Jn 



(4.3) 



Since u is an extremal we have that o'(0) = and we may compute the weak 
variation to be 



(4.4) 



For the inner variation we vary the independent variables, to that end let i : — > IR 
be another function dehned such that 

i(e) = ^[u{x)] (4.5) 

where x = x + e^. Again we have that i'{0) = 0, we then look to compute 

z'(0) = lim i(^[u(5)] - ^[u{x)]) . (4.6) 



We have that 



ix)]^f L{x,uix),^urx)) d~x 



(4.7) 
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where V denotes the spatial gradient operator with respect to x. After a change 
of variable and using the invariance criterion {4.1) we see that 

/ \u{x)\ = ^ \u{x)\ +^Jj'^L -{g^y^^ div(0 + L div(l) + 0(e2) (4.8) 

Subsituting {4.8) into (4.6^ and taking the sum of the inner and outer variations 
concludes the proof. □ 

4.5. Theorem (strong conservation law weak conservation law). Let the vari- 



ational problem (2.9 1 have a variational symmetry in the sense of Definition 3.4 
and that the minimiser to the variational problem is smooth u e C^(f2), then (4.2) 
holds. 



Proof Recall from (3.13) we have that 



(4,10) 



Using Definition 3.3 of the characteristics, Q we see 
which after some simplification yields 

Integrating over the domain concludes the proof. □ 
Now we have developed the framework sufficiently to state our main result in 
this section. Here we are concerned with broken extremals, that is, functions whose 
derivatives have finitely many jump discontinuities. 



4.6. Definition (broken extremal). An extremal, u G C"(r2), to the problem (2.9) 
is said to be a broken extremal if it is piecewise C^(17) over the domain fl with 
bounded derivative. That is, can be decomposed into finitely many open subsets, 
{^i}iLi such that 

(1) the subsets make up the entire domain, i.e., ^ — {Ji ^i, 

(2) they are non-overlapping, i.e., f2,; n fij = and 

(3) the solution is smooth over each of the subsets, i.e., u £ C'^{fli) O W^(r2). 

4.7. Definition (skeleton and jumps). We define 

^ ■.= \J{xedn,} (4.12) 

i 

to be the skeleton of the decomposition. Let rii be the outward pointing normal to 
r2i, we then define jumps of scalars and vector valued functions as 

N := v\nini + u|n2«2 (4.13) 
lvl:={v\n,yn,+iv\n,yn2, (4.14) 

respectively. 



4.8. Definition (piecewise variational symmetry). Let u be a broken extremal to 
the variational problem (2.9) and let {fii}i=i be the decomposed domain of u. Then 
(3.1 ) is a piecewise variational symmetry if it is a variational symmetry over fi^ for 
each i = 1, ... TV. 
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4.9. Theorem (conserved quantities for broken extremals of the variational prob- 
lem). Suppose u is a broken extremal to [2.9) and the problem has a piecewise 
variational symmetry. Then 

N 



- E i + ^ div ^ - ( ^) V. div + / 



dL 



diVu) 



(4.15) 



Proof The proof of this result follows the same lines as the proof of Theorem |4.3[ 
Note that since u is piecewise smooth the Euler Lagrange equations vanish over 
each subdomain fli. Hence we are able to invoke Theorem |4.3| upon each fli. To 
conclude we note that 



/ dL 



V0 



E 

i 

E 



div 



div 



dL 



dL 

a(vw) 



/ dL 

dL 



(4.16) 



9(Vu) 



The other jump term follows similarly by splitting the inner variation from the 
proof of Theorem 4.3 into subdomains. □ 



5. Finite element conservation laws 

5.1. Discretisation. In this section we calculate the discrete counterpart to The- 
orem 4.9 in the finite element context. To that end, let ^ he a, conforming trian- 



gulation of il, namely, =^ is a finite family of sets such that 

(1) K ^ 3" implies K is an open simplex (segment for d—\, triangle for d — 2, 
tetrahedron for d = 3), 

(2) for any K, J & 3 we have that K O J is a full subsimplex (i.e., it is either 
0, a vertex, an edge, a face, or the whole of K and J) of both K and J and 

(3) ijKe^K = n. 

We let be the skeleton (set of common interfaces) of the triangulation and 
say e G ^ if e is on the interior of Q and e € dQ if e lies on the boundary dfl. 
The shape regularity of 3 is defined as 

PK 



inf , 



(5.1) 



where pk is the radius of the largest ball contained inside K and hj,^ is the diameter 
of K. We use the convention where /i : 51 — > IR denotes the meshsize function of 
i.e., 

h{x) -.^ maxhK, (5-2) 
where hn is the diameter of an element K. We introduce the finite element space 



V 



e H^fl) : $|k e P'^VK e ^} n C°(17) n Hj(l7), 



(5.3) 



where P*^' denotes the linear space of polynomials in d variables of degree no higher 
than a positive integer k. We consider fc > 1 to be fixed and denote by N :— dimV. 



The Galerkin approximation to the variational problem (2.9) is to seek ?7 G V C 
Ho(fi) such that 

^im^l^lJiV] (5.4) 



The finite element scheme defined by (5.4) is guaranteed to be well posed under 



some assumptions on L that allow us to invoke the Lax-Milgram Theorem or the 
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more generally applicable inf-sup condition |EG04j . Henceforth we now assume the 
continuous minimisation problem admits a unique solution. 

We may now proceed in deriving a finite element Noether type conservation law. 
As already seen, the conservation law arises after taking inner and outer variations 
of the variational problem. The outer variation can be characterised by the following 
Lemma. 

5.2. Lemma (discrete Euler-Lagrange equations). The discrete Euler-Lagrange 
equations associated to the variational minimisation problem ( 5.4 ) are to seek [/ S V 
such that 







div 



dL 



K 



dL 



d{VU) dU 



V- 



dL 



diVU) 



(5.5) 



Proof Define the real valued function which we call the outer variation operator 



o{e) 



[ L{x,U + eV,\/U + e\/V) 

JQ. 



(5.6) 



where y e V C HJ is the discrete variation. Since [/ G V is the discrete minimiser 
of the energy functional we certainly have that o'(0) = and we may explicitly 
compute this quantity, the first variation. 

dL 



o\e) 



(., U + eV, VU + eVV) )\v + ^^(^. ^ + ^ W + ^VV) ^ 

(5.7) 

Note that since VC/ is not continuous over the skeleton of the triangulation ^ 
we have that gj^jy is, in general, not well defined. But V[/ is smooth over the 
interior of each element. We thus split the integral into elementwise contributions 
and integrate by parts elementwise. For brevity we note that the Lagrangian L = 
L{x, U,VU) and drop the dependency. 



o'(0) 



dL 



dL 



dU 
dL 



V 



IrJ^-'"'''' Oiyu) ' du 



V 



(5.8) 



where div^' denotes an elementwise divergence. We now use the identity 



E 



dL 



diyu) 



and hence 
= 

as required. 



o'(0) 



div 



dL 



K 



d{VU) duJ 



{{V] 



dL 



d{VU) 



V. 



(5.9) 

(5.10) 
□ 



5.3. Example (discrete Laplace's problem). For example the discrete Euler-Lagrange 
equations associated to Laplace's problem (Example 3.9) are to find [/ e V such 
that 

' {AKU + f)v+ f iwuuv}=o yvev, (5.11) 

where A/^ is an elementwise Laplacian. 



Note that if J7 is a piecewise linear function, the first term of (5.11) is zero 



Hence the discrete Laplacian can be completely characterised in terms of the jump 
of the gradient of U over the internal skeleton. 
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5.4. Definition (L2(ri) projection operator). We define Py : L2(ri) — > V such tliat 
for each w E L2(0) we have 

(Pyw,V) = {w,V) VFeV. (5.12) 

5.5. Theorem (conserved quantities over C°(n)-finite element spaces). Let u he 
the unique weak extrema to the minimisation problem (2.9) and U be its finite 
element approximation. Suppose that this problem satisfies a weaic variational 
symmetry. Then the finite element solution satisfies the following 



=J^[U] 



div 



K 



( dL \ dL 



(Vk^)^I + idiv| 
Pv<^l 





■ dL 1 


L 


a(vf/)J 



Mia 

(5.13) 
(5.14) 



Proof Recall we have an energy functional of the form 

^[u] / L{x,u,Vu) . 
Jn 

We have that a finite element minimiser of this energy functional is continuous over 
the domain f2 but its derivative is not. 

Using the discrete Euler-Lagrange equations from Lemma |5 . 2| testing with we 
need only calculate the inner variation. 

The inner variation can be regarded as a change of variables on the independent 



variable jGH961 §3.3]. In a similar calculation to that of the Proof of Theorem 4.3 



we let i{e) = ^ \U(x)\ with x = x + e$,. We again split the integral into subdomains 
to obtain 



0^lun-^{J[Uix)]~ ^[Uix)]) 



lim - 



L{x,U(x),VUix) dx 
+ f LCtik) 

JdK ^ 

/^(V.L)^| + Ldiv^ + (^) Wdiv^ 



L{x,U(x),\IU{x) da;) 



(5.15) 



+ 



Mia, 



as required. 



□ 



5.6. Applications to the p— Laplacian. In this section we give a numerical ver- 
ification to Theorem |5.5| for a simple test problem, that of the p-Laplacian 

-divflVuP"^ Vw) 



/, (5.16) 
where we will restrict p G (l,oo) for degeneracy issues. The p-Laplacian is the 

o 

Euler-Lagrange equation of the following minimisation problem: Find u e Wp(il) 
such that 

AM<AM V«eWi(0), (5.17) 
with the (parametrised) action functional given by 



J a V 



(5.18) 
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Note that for p = 2 this problem coincides with the standard Laplace's problem 
(see Example 3.9 1. For general p it is well known that the problem is uniquely 
solvable. 

The discrete weak formulation associated to the minimisation problem (5.171 is 
to find J7 e V such that 



/ |v;7|^ ^(vc/)'^VF= ( fv vv^ev. 

Jn Jn 
In this test we choose / such that 



u — sm TT a; 



(5.19) 



(5.20) 



solves the p-Laplace equation (5.16). We have that / can be written as f — f{\x\) 
and hence the Lagrangian 



L(x,v,\/v) = - iVul' 
P 



fv 



(5.21) 



is invariant under SO{d) group actions. 

We fix d = 2 and take 5^ to be a structured triangulation of il, the unit circle, 
as given in Figure [T] 



Figure 1. An example of the triangulation and the finite 
element approximation to m = sin ^tt |a;|^), the solution of the p- 
Laplacian. 



.:„--..H;.:i.":i.j.-*--v-*^-'^ 




(a) An example of the triangulation, here (b) The piecewise linear finite element approxi- 
dim V = 12564. mation of the solution of the 3-Laplacian. 



It is well known |BL941 c.f.] that the finite element approximation (5.19) is 
well posed and has optimal convergence properties. In Tables [T]-[3] we show errors, 
convergence rates and the values of the finite element Noether quantity as written 
in Theorem |5.5| for various cases of p. The tables also study the experimental order 
of convergence of the numerical approximation which we now define. 

5.7. Definition (experimental order of convergence). Given two sequences a{i) and 
h{i) \ 0, i — I, . . . „ we define experimental order of convergence (EOC) to be the 
local slope of the loga(z) vs. logh{i) curve, i.e., 



(5.22) 
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Table 1. In this test we computationally study the behaviour 
of the finite element conserved quantity, ^[U], given in Theorem 
5.5| To that end we consider the p-Laplacian with p = 3. We 
fix / such that u is known and is given by (5.201. We compute 
the piecewise linear (fc : 
(p9|. 



1) finite element approximation given by 
To do this we formulate (5.19) as a system of nonlinear 



equations, the solution to this is then approximated by a Newton 
method with tolerance set at 10^^°. At each Newton step the 
solution to the linear system of equations is approximated using 
a stabilised conjugate gradient iterative solver with an algebraic 
multigrid preconditioner, also set at a tolerance of 10"^". Note 
that y^[U] is below the tolerance of the solvers when the mesh is 
sufficiently resolved. 



dimV 


Ir ^ IIlp(o) 


EOC 


II iii,p 


EOC 


,yr[U] 


13 


1.12725161 


0.000 


3.63060492 


0.000 


2.353673e- 


14 


41 


0.70442091 


0.678 


2.98329491 


0.283 


1.776357e- 


14 


145 


0.15390246 


2.194 


1.56695787 


0.929 


3.108624e- 


15 


585 


0.03539738 


2.120 


0.73020113 


1.102 


1.526557e- 


13 


2805 


0.00618342 


2.517 


0.31109110 


1.231 


2.456799e- 


12 


14293 


0.00110803 


2.480 


0.13244566 


1.232 


2.600439e- 


12 


73401 


0.00022205 


2.319 


0.05750230 


1.204 


6.134856e- 


13 


384769 


0.00004666 


2.250 


0.02645758 


1.120 


1.581680e- 


13 



Table 2. This is the same problem as in Table [T] with the excep- 
tion that we take p — A. 



dimV 


\\u - U^W 


EOC 


II iii,p 


EOC 




13 


1.50645187 


0.000 


3.89649998 


0.000 


4.002354e-14 


41 


0.95656811 


0.655 


3.32660211 


0.228 


1.776357e-14 


145 


0.18567964 


2.365 


1.69441051 


0.973 


3.730349e-14 


585 


0.04346809 


2.095 


0.77624620 


1.126 


1.314726e-12 


2805 


0.00787713 


2.464 


0.33055071 


1.232 


2.027223e-ll 


14293 


0.00141755 


2.474 


0.13947531 


1.245 


2.106114e-ll 


73401 


0.00028203 


2.329 


0.06028734 


1.210 


4.074260e-12 


384769 


0.00005934 


2.249 


0.02756628 


1.129 


1.291921e-12 



Table 3. This is the same problem as in Tabic [T] with the excep- 
tion that we take p = 5. 



dimV 


11" ^ llLp(O) 


EOC 


II lll,p 


EOC 




13 


1.84301273 


0.000 


4.16722611 


0.000 


7.812362e-13 


41 


1.14195690 


0.691 


3.69726521 


0.173 


8.526513e-14 


145 


0.21243009 


2.426 


1.85724530 


0.993 


2.984279e-13 


585 


0.05034349 


2.077 


0.83769960 


1.149 


1.120704e-ll 


2805 


0.00932480 


2.433 


0.35992898 


1.219 


1.558869e-10 


14293 


0.00170255 


2.453 


0.15093728 


1.254 


1.587779e-10 


73401 


0.00033701 


2.337 


0.06531226 


1.209 


2.558158e-ll 


384769 


0.00007080 


2.251 


0.02955021 


1.144 


1.013448e-ll 
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5.8. Mimetic methods weakly enforce discrete conservation laws which 
are derived from trivial Lie group actions. The mimetic finite element frame- 
work consists of reformulating the Euler-Lagrange equations as a system of first 
order PDEs. Consider our prototypical example for illustrative purposes. Poisson's 
problem, 

Au = 0, (5.23) 
is the Euler-Lagrange equation of the minimisation problem 

^[u]= / ;^ |Vu|^ ^ min. (5.24) 

It can be written in mixed form by introducing an auxiliary variable p to represent 
the gradient and rewriting Poisson's problem to seek(w,p) such that 

divp = (5.25) 

p = Vu. (5.26) 

These are then the Euler-Lagrange equations of the saddle point problem 

je[u,p] := / l\p\^ +uidivp). (5.27) 
Jn ^ 

The correct function space setting is to seek u G L2(fi) and p G H'^'^(O) 
{^f : div^ G L2(r2)}. A conformal approximation of this problem can be sought 
using the Raviart-Thomas and piecewise constant finite element pair ^RTTTj, for 
example. A sufficient condition for the construction of a conformal finite element 
space ofH'^'^(rj) is that the jumps of the discrete functions vanish over the skeleton 
of the domain. Taking a generic P G IRTq we see for our model problem 



= {div p\(i>) = - P'^V4> 

Jn 



riKq 



dK 



(5.28) 



E / divifP0- E -P"^' 



K^sr-^^ eean 
This gives a "weak" Stokes theorem. 

Recall Remark |3.10| concerned itself with the trivial Lie group action of trans- 
lation in the dependent variable. For our model problem we have that Vu is a 
conservation law. The mimetic scheme weakly enforces this conservation law. 

6. Conservative properties of Lagrangian FEs for strong solutions 

In this section we present results concerning the approximability of the strong 
continuous conservation laws arising from Theorem |3.6| We examine numerically 
the behaviour of the Lagrangian finite element method. In this sense we wish to 
measure the quantity div('^[t/]) and evaluate how far it deviates from zero. For 
clarity of exposition, we will assume henceforth that the continuous minimisation 
problem takes the form: Find m G Hq such that 

J\u\^ inf ^H. (6.1) 

v&A\{n) 
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6.1. Theorem (Bound on the finite element approximation of Noethe r's l aws.). 
Let u e H^(f7) n Hj(il) be a strong extrema to the variational problem {2.9) (and 
hence a strong solution to the Euler-Lagrange equations { 2.1^ ). Suppose we have 
that Theorem 3.6 holds under a variational symmetry group with inhnitesimals ^ 
and (j). 

In addition assume that we have 



dL{x, u, Vii) 



(6.2) 
(6.3) 



Then if U G V is the finite element approximation to u there exists a constant C 
such that 



|div'^[J7]||jj-i(o) < Cy\\L{x,U,\/U) - L{x,u,\/u)\\ + ||VC/- Vi 

dL{x,U,VU) dL{x,u,Vu) 
d{S7U) d{Wu) 

^{x,U)-(t>ix,u)\\ + Uix,U)-ax,u) 



(6.4) 



Proof We begin by noting that since m is a strong extremal, Theorem 3.6 holds 
and we have that div'^[u] = 0. Hence 



l|div^[[/]|lH-i(a) = WdivmU] - ^M)llH-(n) • 
Now, we may use the fact that for a generic G HJ 

<mu]-nu]\\\\vv\\- 



(6.5) 



(6.6) 



Since (p was generic we may divide through by ||V(/3|| and take the supremum over 
if then by the definition of the H^^(ri) norm 



||div^[[/]|lH-i(o)<ll'^[c/]-'^MII 



(6.7) 



By the definition of the Noether quantity from (3.10) we have that 

dL{x,U,VU) dL{x,u,Vu 



|div^[C/]|lH- 



< 



d{WU) 



{x, u)- 



(|(a;,C/))^Vt/ — ^ — - -{i{x,u)y\/u ^ ' ' 



d{VU) 

+ \\L{x, U, VU)^{x, U) - L{x, u, Vu)^{x, u 



a(Vu) 



(6.8) 



where for clarity we have written the dependencies explicitly. Now for each of 
the J^i we add and subtract appropriate quantities and make use of the triangle 
inequality. We thus have the following bounds: 



{x,U) 



dL{x,U,VU) dL{x,u,Vu) 



d{VU) 
i(j}ix,U) ~ (t)ix,u)) 



9(Vu) 
dL{x, u, Vu) 



(6.9) 
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< 



{^ix,u)-^ix,U)){Vu 



J dL{x, u, Vm) 



i^{x,U)y{Vu-VU) 

{^{x,u)yvu 



a(Vu) 

dL{x, u, Vw) 



and 



dL{x,U,VU) dL{x,u,Vu) 

d{vu) a(Vw) 



^3 < \\{L{x,U,VU) - L{x,u,Vu))i{x,U)\\ 
+ \\L{x,u,Vu){S,{x,U) ^ i{x,u))\\ . 



(6.10) 



(6.11) 



Under assumptions (6.2)- m.M and using the fact that and ^ are infinitesimals 



of a smooth Lie group action we have that 

dL{x,U,VU) dL{x,u.yu) 



=^1 < C 



d{VU) d{Vu) 
■^2 < C{ U{x, u) - ^{x, U)\\ + \\Vu - VC/| 



\-mx,U)-(bix,u)\\) (6.12) 
dL{x,U,VU) dLix,u,Vu) 



a(v[/) a(vu) 

^3 < C{\\L{x, U, \/U) - L{x, u, Vu)\\ + U{x, u) - ^{x, U)\\) . 
Taking the sum of the J^i gives the desired result. 



(6.13) 
(6.14) 

□ 



6.2. Corollary. Let the conditions of Theorem \ 6.1\ hold under the same variational 
symmetry group with inhnitcsimals ^ and 4>. In addition assume that the La- 
grangian is sufRciently smooth such that both L and Qm^\ are (locally) Lipschitz 



with respect to the second and third variable then the bound {6.4) can be simplihed 
to 



l|div<^[C/]|lH-i(a) <C\^\\U-u\\ + \\VU - Vu\\ 

+ Uix, U) - c^{x, u)\\ + \\£,[x, U) - £,{x, u)\\ 



(6.15) 



6.3. Remark (relating to the p-Laplacian) . We may relate Theorem 6.1 and Corol- 
lary |6]2] to the p-Laplacian studied in |5.6[ We were considering that L was invari- 
ant under rotations in the independent variable. In that case we have that = 
and that 



hence we may infer that 



Idiv^[[/]||H- 



U{x,U)~£,{x,u)\\^Q 



< C \\U-u\\ + live/- Vull 



(6.16) 



(6.17) 



The leading term here is ||VC/ — Vw|| which is well known to be = 0{h'^). 

6.4. Lemma (a computable upper bound for the conse rvat ion law). Let u G H^(f2)n 
Hg be a strong extrema to the variational problem {2.9). Suppose that Theorem 
\3.(^ holds and that U is the finite element approximation to u. Then there exists a 
constant C dependent on the shape regularity of 3" such that 



l|diV^[C/]|iH- 



n-Hn) < mJ) C-f E l|diV^[C/]|lL,(K) + E li I'^If^H llL.(e) 



(6.18) 
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Proof The proof is a standard aposteriori argument where the quantity of interest 
is spht into regular and singular parts. Recall for the model elliptic problem "^[C/] = 
'^{x, U, VU) and hence div"^[;7] ^ L2(rJ). Let (• | •) denote the R'^n) ~ HJ duality 
pairing then for any </3 e Hq it holds that 



{div'^[U]\ip) ^ J2 ^ [ '^l^Y'^ip 

Ke.9' 



r r (6.19) 

^^^jk J ok 



Note that for generic p g [L2(il)]'' and v e H^(0) that 

y2 / p^riKV = / [p1 w + / p^riKV. (6.20) 

K<z,'^-'dK Js Jan 

Hence 

(div<r[c/]i^)- ^ / dw^^[u]^-j2 [ mmif- E hpVn^- (6.21) 

Applying a Cauchy-Schwarz inequality followed by a Poincare inequality together 
with a trace inequality yields 

(div<r[C/]|^) <C l|div^[f/]llL.mll¥'llL.W+EliI"^[^]lllL.(e)ll^llL.(e) 



<C||V^|| ^ ||div^[[/]|lL,(;,) 



(6.22) 

Noting iy9 was a generic function in Hq, dividing through by HV^sjl and taking the 
supremum over ip yields the desired result. □ 

6.5. Remark. It is worth noting that since the quantity div'^[C/] is not orthogonal 



to V there are no powers of h appearing in the estimate given in Lemma 6.4 What 
is interesting is that the estimate still converges to zero at the same rate as that of 
the residual of the problem. 

6.6. Numerical experiments. In Figure[2]we show numerically that the estimate 
given in Lemma [6.4| converges at the same rate as the apriori bound given in Remark 
|6.3[ All numerical experiments are conducted for simplicity on the 2-Laplacian 
taking = [— 1, 1]^ which is discretised using an unstructured triangulation. 

6.7. Adaptive methods. We conclude our numerical experiments in Figure [3] by 
using the observed convergence of the computable estimator E{U, f) to construct 
an adaptive scheme aimed at minimising E{U,f) (and hence ||div'^[t/]||_j^). The 
adaptive algorithm we make use of is of standard type (SOLVE — ^ ESTIMATE 
— )■ MARK — REFINE |SS05[ c.f.]) utilising the maximum strategy marking and 
newest vertex bisection refinement. 

7. Conclusions and outlook 

In this work we have given a consise statement of Noether's first theorem, applied 
to a set of model problems. 

We have proved a Noether type theorem for a specific class of weak extremum 
to a model variational problem, that is, those that are W^(r2) with finitely many 
jump discontinuities. 
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Figure 2. In this experiment we consider the 2-Laplacian. We 
fix / such that u = exp(^— 10 \xfj. We solve the discrete problem 
on concurrently refined meshes and compute the L2(r^ )-er ror, the 
Hp-error and the computable estimate given in Lemma i 
that in each of the examples E{U, f) converges like 0(7i 



6.4 



Notice 
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(a) Initial mesh 



(b) In this example we choose k - 
dim V = 295681, £;([/,/) 0.08. 
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"^EOC = 3.99 



(c) In this example we choose k = 2. Note when (d) In this example we choose fc = 3. Note when 
dimV = 295681, £{(7,/) as 0.0003. dimV = 295681, £((7, /) as 0.000001. 



We have in addition proved an equivalent discrete theorem for the finite element 
approximation of the problem. We write the exact conserved quantity for the 
discrete scheme in the same spirit as |,HM04]. 

We have demonstrated that the Lagrangian finite elements enjoy the property of 
asympotically conserving the strong Noether conservation laws when approximating 
strong solutions of certain classes of variational problem and Lie group action. 

We have studied the exact discrete conserved quantities numerically and in ad- 
dition also constructed a geometric based adaptive scheme to conserve the approx- 
imate continuous quantities up to user specified tolerance. 

In a forthcoming sibling publication [MP12j we will examine the structure of the 
conservation laws for less regular finite element spaces, the discontinuous Galerkin 
methods. In this case the arguments are generalisable to higher order Lagrangians. 
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